Fractional Energy Loss and Centrality Scaling 
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The phenomenon of centrality scaling in the high-pr spectra of 7r° produced in Au-Au collisions 
at ^/i = 200 GeV is examined in the framework of relating fractional energy loss to fractional 
centrality increase. A new scaling behavior is found where the scaling variable is given a power-law 
dependence on JV par t. The exponent 7 specifies the fractional proportionality relationship between 
energy loss and centrality, and is a phenomenologically determined number that characterizes the 
nuclear suppression effect. The implication on the parton energy loss in the context of recombination 
is discussed. 
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PACS numbers: 25.75.Dw 

The production of hadrons in Au-Au collisions at the 
Relativistic Heavy-Ion Collider is found at high pr to 
depend sensitively on centrality 0. At y/s = 200 GeV 
the scaled inclusive cross section of pions at pr ~ 3-4 
GeV/c decreases by a factor of 4-5 when the centrality 
is varied from the most peripheral to the most central 
@,0|. In a previous paper we reported the finding of a 
universal function <f>(z) that can describe all the inclusive 
cross sections at all centralities. Such a scaling behavior 
is achieved by use of a scaling variable z that combines 
Pt with iVpart, the number of participants. In this paper 
we investigate the origin of that scaling. In particular, 
we consider the nature of energy loss that can give rise 
to such a behavior. 

Although the scaling behavior can be extended to in- 
clude energy dependence also p|, we restrict our con- 
sideration here to centrality dependence only, and em- 
phasize the phenomenological implication of the data at 
s/s = 200 GeV, measured by PHENIX 0. We shall not 
discuss energy loss at the parton level except near the 
end, since perturbative QCD is not reliable for pt < 6 
GeV/c. Indeed, soft partons have been found to be im- 
portant in the hadronization process through recombina- 
tion [1, II, @ ■ We shall stay mainly at the hadronic 
level that is phenomenological and uncontroversial, and 
consider the nuclear suppression effect on the observed 
hadrons. 

Let us recall the scaling behavior found in With 
the definition of the variable 



z = p T /K{N) 
we find that the function $(z), 

1 dN v 



= A{N)K\N) 



2irpT drjdpx ' 



(1) 



(2) 



exhibits scaling behavior. We have used the notation 
N = -/Vp ar t, for brevity, and 

K(N) = 1.226 — 6.36 x 10~ 4 A^ , (3) 
A(N) = 530N c (N)-°- 9 , N C (N) = 0.44JV 1 - 33 , (4) 



where both K(N) and A(N) are normalized to 1 at N = 
350. At all pt and N, the data of dN^/pTdrjdpT at 
midrapidity collapse to one universal curve ^(z), which 
is parametrized in Q by 



$(z) = 1200 (z l + 2)~ 4 ' 8 (1 + 25e" 4 - t>z ) 



(5) 



It is clear that what gives rise to the scaling behavior 
must be related to an universal property in the medium 
effect on the production of pions. Since it is not possible 
to determine experimentally the degradation of parton 
momentum as the medium size is increased, and since 
whether hadronization is by means of fragmentation or 
recombination is still controversial, we choose to stay 
at the hadronic level and examine energy loss. Since 
the produced pions do not themselves traverse the dense 
medium, energy loss here does not refer to the evolution- 
ary process of a pion, as one can for a parton. Instead, it 
refers to the shift of the pion distribution, as the medium 
size quantified by N is increased. 

Let us now consider the implications of a scaling func- 
tion $>(z). Let z be defined as 



z = xJ(N) , 



(6) 



where a; is a dimensionless momentum variable identi- 
fied as x — pt/po with po chosen at po = 1 GeV/c so 
that x is numerically the same as px- We now examine 
the consequences of writing <fr(z) in terms of x and N 
explicitly 



= F(x,N) 



(7) 



In pQCD, such as in Ref. [13], one compares a distribu- 
tion in medium with one in vacuum in order to emphasize 
the medium effect. We prefer, however, to stay away from 
the pp collision case, since we want to consider incremen- 
tal changes of the medium size. From our perspective of 
dealing only with the observables, it is very natural to 
ask the e-S type question. That is, given a medium that 
is not too small, in which pions are produced, if its size 
is increased by an e amount, what is the corresponding 
downward shift 5 in momentum in order to maintain the 
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same probability of producing the pions? In terms of x 
and N, the proposition can essentially be stated as 



F(x,N) = F(x-S,N + e) 



(8) 



The +e and —S relationship is a consequence of the sup- 
pression effect. 

For infinitesimal e = SN and 8 = Sx we can expand the 
right-hand side of Eq. JSJl and keep only the first order 
terms, getting 



Sx _ dF/dN 
SN ~ dF/dx ' 

If F(x, TV) satisfies Eqs. © and J7J), then we have 

Sx x dJ 
SN ~ JdN ' 



(9) 



(10) 



Since Sx is proportional to x, it is more sensible to con- 
sider the fractional energy loss 



Sx 
x 



din J SN 
din TV TV 



(11) 



If the fractional energy loss is proportional to the frac- 
tional change of centrality, a notion that seems extremely 
reasonable, i.e., 



Sx _ SN 
~x ~ 7 ~N 



(12) 



where 7 characterizes the suppression effect, then it is 
necessary that 



J(N) = 



N 



(13) 



for some normalization Nq. Clearly, Eq. (|12() does not 
make sense for N very small, such as N = 2, since e 
cannot be made infinitesimal compared to 2. 

The power-law behavior in Eq. (|13[) is a necessary con- 
sequence of scaling and fractional proportionality, Eq. 
(jT2|l . Comparing J (TV) with 1/K(N) in Eq. ©, one 
finds that Eq. (|13f) differs enough from our first scal- 
ing parametrization to cast some doubt on whether Eq. 
(|13fl is sufficient for scaling. However, the fractional pro- 
portionality relationship is so compelling that we have 
been motivated to reexamine the data, especially since 
the original preliminary data have by now been finalized. 
It should be noted that Eq. ljT3|) is obtained without re- 
lying on any specific form for $(2;); it depends only on 
the structure of the scaling variable expressed in Eq. JSJl . 
Thus one expects Eqs. ljr2")l and (fH5|) to be very general 
properties of centrality scaling. 

Using the data on tt° from PHENIX tabulated on the 
web 111 we have assembled the spectra for all measured 
centralities, and made appropriate horizontal and vertical 
shifts in the log-log plot to obtain a universal behavior. 
The result is shown in Fig. 1. The ir + data for the most 
central collisions are used to supplement ir° in the low-px 



10 



N 

<6> 



10" 



° 7t + 0-5 % 

* 71° 0-10% 

► ji° 1 0-20% 
< jt° 20-30% 

* Ji° 30-40% 

* tc° 40-50% 

* jc° 50-60% 
« 7t° 60-70% 
- jc° 70-80% 

* 7i° 80-92% 



10 



10 



FIG. 1: Scaling distribution $(z) showing the coalescence of 
all data points of 9 centrality bins for tv° and the most central 
collisions for ir + production in Au+Au collisions at y/s = 200 
GeV as measured by PHENIX [Hill. The solid Hue is a fit 
parametrized by Eq. I)15|l . 



region 12]. The horizontal shift determines the scaling 
factor J(N) shown in Fig. 2(a). The vertical shift deter- 
mines the normalization factor B(N) shown in Fig. 2(b). 
The resultant scaling distribution $(z) is now related to 
the measured inclusive distribution by 



B(N) 



Po 



dN„ 



J 2 (N) 2-kpt drjdpr 



(14) 



Evidently, the scaling behavior exhibited in Fig. 1 is very 
good. The solid line is a fit using the formula 



$(z) = 150 (z 2 + 1.05)" 



(1 + 7e- 4z ) . 



(15) 



There are many more points included in Fig. 1 than those 
of the preliminary data used in Ref. Q|. 

The behavior of J(N) in the log-log plot in Fig. 2 can 
be fitted by a straight line according to Eq. $T3$ with 



7 = 0.077 



(16) 



The region TV < 10 is not considered. The normalization 
point is chosen to be TVo = 325, which corresponds to 
the most central 0-10% collisions 0- The normalization 
factor B(N) can be fitted by two straight lines as shown 
by the solid lines. The implication of this result will be 
discussed in the following. 

Although Eqs. i|13|) an d (|ltjfl do not differ too much 
from the inverse of K(N) given in Eq. @, we have in- 
vestigated the source of the difference. Our conclusion 
is that in we read the preliminary data from a figure 
given in a conference talk |13( and estimated the central 
points, whereas here we use the finalized data in tabu- 
lated form, which differ slightly from the original. Thus 
our present result is more reliable. Moreover, since in our 
new description the medium suppression effect is charac- 
terized by one and only one parameter 7, which plays 



3 



1.1 
1 

— 0.9 
~~> 0.8 

0.7h 

2 



10 



CO 



10' 



10" 



(b) 



(x) is related to the inclusive cross section as 
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FIG. 2: (a) The scaling factor J{N). Note that the vertical 
scale is logarithmic. The solid line is a straightline fit, yielding 
7 = 0.077. (b) The normalization factor B(N). The solid 
lines are fits showing two scaling regions. 



the crucial role of specifying the fractional proportional- 
ity relation in Eq. i|I2|) , the resultant scaling behavior has 
the distinction of being physically motivated. As we have 
shown, it is not so much what the scaling function ^(z) 
is as what the scaling variable z is. Since z quantifies 
the difficulty of producing transverse motion, it can be 
termed transversality that gives a universal description 
of that difficulty at any centrality. The z dependence of 
<£>(z) in Eq. ljT5|> differs little from that of Eq. J5J), except 
at very high z where the new data extend beyond that 
of the preliminary result. 

With the value of 7 in Eq. (|16|) now determined phe- 
nomenologically, we can return to Eq. (|12(l and claim that 
the notion of fractional proportionality has direct support 
by the data. Note that the independence of the fractional 
energy loss Sx/x on a; is a result that differs from that 
of pQCD on the momentum shift of hard partons where 
Apt is proportional to p^ 2 [TH ]. and is consistent with 
the independence of Raa on px • Equation Ijl2(l also 
implies that Sx/SN is proportional to x/N which is very 
different from the assumption that the energy loss per 
unit length traversed by a parton, dE/dL, is a constant. 
Although our phcnomcnological result on the produced 
pions has no direct implication on the evolutionary prop- 
erties of the partons propagating through a medium, one 
should keep these differences in mind, when the measur- 
able consequences of what can be calculated in pQCD 
are inferred. 

Since the method by which we obtained J(N) is by 
means of data fitting, it is desirable to find an alterna- 
tive method that is more direct. Given $(z), one can 
determine the average 

»-/* 0.42. (17) 

Then from Eq. © we have J(N) = 0.42/ (x) (TV), where 



(x) (N) = / dp T PT 



dN„ 



Po / dpj 



dN n 



(18) 



drjdpx I J drjdpT 

Since J(N) is normalized to 1 at N = No, we now have 

J(N) = (x) (N )/ (x) (N) , (19) 

thus eliminating the reference to the scaling variable z. 
The N dependence of (x) (N) can be determined directly 
from the data when the average is calculated at various 
centralities. The large px part of the integration is impor- 
tant, for it is that part of the dN^ / drjdpx that led us to 
the scaling factor J(N) in the first place. The 1ow-pt part 
of the spectra should be accurately parametrized. Now, 
it is a matter of evaluating Eq. HI8(I instead of shifting 
and rescaling in Fig. 1. We recommend that experimen- 
tal groups that have the data, not only of tt° , but also of 
charged hadrons, can use this method to check whether 
J(N) indeed has the form given in Eq. I|I3I) . 

The scaling behavior that we have obtained is for the 
produced 7r°. We have found in that the anoma- 
lously high p/tt ratio at px ~ 3-4 GeV/c can be under- 
stood in terms of recombination without fragmentation, 
when hard partons are allowed to recombine with the 
soft ones. That is possible if the recombination function 
for the pion does not restrict the recombining q and q to 
have roughly the same momentum [l4| . Indeed, in our 
view jet fragmentation is included in recombination be- 
cause a large-pT hard parton initiates a parton shower 
that hadronizes by recombination 0] . The use of frag- 
mentation function is only a phenomenological way of 
parametrizing that process, and does not stand for an 
independent hadronization mechanism. In heavy-ion col- 
lisions the parton shower on the surface has low-px com- 
ponents that mingle with the soft partons with hydrody- 
namical origin. Since separating them would be artificial, 
we treat them on equal footing in the recombination pro- 
cess that can involve hard partons as well. The purpose 
of this discussion is to prepare our way to descend to 
the parton level without fragmentation. The scaling be- 
havior that we have found supports this view, since our 
rescaling procedure is universal and does not separate the 
high-pT from the low-px regions using different transver- 
sality variables. We note that this view differs from that 
taken in Q, in which recombination and fragmentation 
are important in different regions. We also note that 
our application of the recombination mechanism to all 
partons encounters no inconsistency with the two exper- 
imental facts: (a) p/ir ratio is roughly I at px ~ 3-4 
GeV/c, and (b) the jet structure in Au-Au collisions is 
similar to that in pp collisions in the same px region 

With the above discussion we have laid the basis for 
the expectation that the quark distribution F q for all px 
contributes to the determination of the pion distribution 
at all px- The recombination equation derived in Refs. 
[H |(| has the form 

$(z) = J ' dz 1 dz 2 ^F q (z 1 )F q (z 2 )8(z l +z 2 -z) , (20) 
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in which the px variables have all been transformed to 
the scaling variables with 

z l = ^J{N), i = l,2. (21) 
Po 

The S function comes from the recombination function 
and guarantees the conservation of momentum. Note 
that z\ and z 2 are integrated over all values; in partic- 
ular, they are not restricted to the region z/2. Indeed, 
since the quark distributions fall rapidly with Zi it is nec- 
essary for one to be large and the other Zi to be small 
in order to give the highest contribution to $(z) at large 
z. The essential remark we want to make now is that in 
our formalism for hadronization the q and q distributions 
must have their pi T scaled by the same factor as shown 
in Eq. l(2T)l. Since the J(N) in Eq. i(2"TJl is the same as 
that for the pion, we conclude that the fractional energy 
loss for the quarks (and antiquarks) satisfies the same 
proportionality relationship as in Eq. i|12|) . It should im- 
mediately be emphasized that these partons are at the 
end of their evolution (hydrodynamical and/or branch- 
ing in showers) just before recombination. As mentioned 
earlier the energy loss discussed here does not refer to 
the radiative energy loss of a parton traversing a medium 
[l7| . Our point is that the fractional energy loss of the 
partons at hadronization satisfies the same property as 
for the pions. Since in this formalism of hadronization 
we have been able to obtain the correct p/ir ratio [6J, it 
follows that the produced protons should have the same 
property in fractional energy loss also. This prediction 
should be checked experimentally by studying the pro- 
ton spectra and seeing whether centrality scaling can be 
achieved with the same J(N). We expect, however, the 
proton mass effect to break the scaling at low z. 

Our final remark is a speculative one. The universality 
of the single exponent 7 for all centralities (except the 
very peripheral collisions) raises the question whether a 
hot and dense medium is any different from a less dense 
medium in its effect on pion production at high px- If 
not, the high px spectra would not be a fruitful place 
to find the signature of plasma formation. A possible 



escape from that conclusion is the observation that the 
result on B(N) in Fig. 2(b) can be fitted by two straight 
lines, which can be parametrized by 

B(N) = {N/Ni)-? 1 , N < 38, 

= {N/N 2 )~ 02 , N > 38 , (22) 

where (5 X = 0.744 and fa = 1.292 with Ni = 1610 and 
N 2 = 325. It suggests that there are two regimes of N, 
requiring different exponents (3\ and (32 to achieve scaling. 
If this break can be ascertained by more detailed analysis 
at different energies, then it offers a way out of the strict 
universality in which there is no hint of any essential 
difference between hot and cold media, leaving no room 
for any signal for the formation of quark-gluon plasma. 

To summarize, we have found a relationship between 
the centrality scaling behavior of the observed pion spec- 
tra and the fractional energy loss of the pions that is 
proportional to the fractional centrality change. That 
proportionality is specified by an exponent 7, which char- 
acterizes the medium suppression effect. In the recombi- 
nation model the same value of 7 is valid for the fractional 
energy loss of the light quarks just before hadroniza- 
tion. The existence of the universal scaling function of 
transversality suggests that there is no essential differ- 
ence in how the low- and high-p*r hadronization processes 
should be treated. The possibility of a single exponent 7 
that can be directly extracted from the data to summa- 
rize the nuclear suppression effect offers a very succinct 
description of a complicated dynamical process. Univer- 
sal behavior in centrality may be broken by the existence 
of two scaling regions in the normalization factor, thus 
providing the possibility of a threshold for a new distinc- 
tive regime. 
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